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ABSTRACT 

A comparative study of the gain achievement is performed in a waveguide optical amplifier whose active 
layer is constituted by a silica matrix containing silicon nanograins acting as sensitizer of either neodymium 
ions (Nd^+) or erbium ions (Er®+). Due to the large difference between population levels characteristic times 
(ms) and finite-difference time step (10“^^s), the conventional auxiliary differential equation and finite-difference 
time-domain (ADE-FDTD) method is not appropriate to treat such systems. Consequently, a new two loops 
algorithm based on ADE-FDTD method is presented in order to model this waveguide optical amplifier. We 
investigate the steady states regime of both rare earth ions and silicon nanograins levels populations as well as 
the electromagnetic field for different pumping powers ranging from 1 to 10^ mW.mm'^. Furthermore, the three 
dimensional distribution of achievable gain per unit length has been estimated in this pumping range. The Nd^+ 
doped waveguide shows a higher gross gain per unit length at 1064 nm (up to 30 dB.cm'^) than the one with 
Er®+ doped active layer at 1532 nm (up to 2 dB.cm'^). Considering the experimental background losses found on 
those waveguides we demonstrate that a significant positive net gain can only be achieved with the Nd®+ doped 
waveguide. The developed algorithm is stable and applicable to optical gain materials with emitters having a 
wide range of characteristic lifetimes. 

Keywords: ADE-FDTD, waveguide, rare earth. Si nanostructures. Gain, Computational methods, optical am¬ 
plifier, silicon photonics 


1. INTRODUCTION 

For many years, rare earth ions have been used in silica-based optical amplifiers. In these systems, the low 
gain value requires to employ significant length (10 to 15 m) of doped fiber to achieve a workable power operation. 
In more compact systems such as erbium-doped waveguide amplifiers (EDWA) a higher gain has to be reached 
in order to shorten the operating length of the amplifier.lll Mostly trivalent erbium ions (Er^+) have been studied 
due to their emission wavelength at 1532 nm, which is adapted to the telecommunications window in optical 
fibers.^ However, there are three major gain limiting factors for the erbium ions : the up-conversion mechanism, 
the excited state absorption and the re-absorption of the signal from the fundamental level. This last drawback 
is characteristic of a three levels system. More recently, neodymium ions (Nd®+) has been proposed instead of 
erbium ion since its four levels emission scheme prevent signal re-absorption from the fundamental level and 
makes it more suitable for achieving higher gain. 

One main drawback of seeking gain with rare earth ions is their low absorption cross section ((Tabs)- However 
this can be overcome by the use of sensitizers that are characterized by a larger absorption cross section and 
an efficient energy transfer to RE^+ ions. Several sensitizers of RE have been proposed in literature, Polman et 
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afil show the sensitization of Er®+ ions by different kinds of sensitizers such as ytterbium ions, metal ions and 
silicon nanograins (Si-ng). Several studieJ3El have pointed out that silicon nanograins are efficient sensitizers of 
Er®+ ions and can increase its effective absorption cross section by a factor up to 10^. MacDonald et aP showed 
a likewise efficient sensitization of Nd®+ ions by those Si-ng sensitizers. 

We aim to model the propagation of an electromagnetic field into a waveguide with a layer containing 
absorbing and emitting centers as for example rare earth ions and silicon nanograins. More particularly, we want 
to determine the system characteristics as fields, level populations, and gain in a steady state regime as a function 
of initial parameters such as concentration of emitting centers, geometry, pumping configuration and pump and 
signal powers. Moreover, to model those steady states regimes in a waveguide with a layer containing absorbing 
and emitting centers, we must take into account the time evolution of electromagnetic field and electronic levels 
populations of these centers. In a waveguide containing silicon nanograins and rare earth ions, the typical lifetime 
of the electronic levels is about some ms, whereas the characteristic period of the electromagnetic field is of the 
order of fs. The steady state of the system cannot be reached with such a difference between absorbing and 
emitting centers lifetimes and electromagnetic field period by means of classical ADE-EDTD. It would require 
prohibitively long computation times with about 10^® iterations. Consequently, we developed a two loops new 
algorithm based on ADE-EDTD method which allows to compute the electromagnetic field distribution, the 
electronic levels populations, and the gain in the waveguide in the steady state regime00 In this paper we 
present briefly the two loops ADE-EDTD algorithm of calculation and the comparative study of waveguides 
with an active layer containing Si-ng and doped either with erbium or neodymium ions. 

2. RARE EARTH DOPED SILICON BASED WAVEGUIDE 

The system describe hereafter is an optical amplifier strip loaded waveguide, composed of three layers as 
presented on figure Bottom and top cladding layers of pure silica of a thickness equal to 3.5 /rm and 0.4 /im 
respectively ensure vertical guiding of modes, whereas a top silica strip stacked on the active layer ensures lateral 
confinement of modes. The static refractive index (i.e. refractive index which remains constant with wavelength) 
of the active layer is equal to 1.5 and the one of the strip and bottom cladding layers to 1.448 to ensure the 
guiding conditions. We inject at one end of the waveguide pump and signal guide modes in a co-propagation 
scheme as presented on figure 1(a). The active layer of 2 fj,ia thick is constituted of non-stoichiometric Silicon 
oxide (SiOx) containing silicon nanograins and rare earth ions (RE®+). 

The active layers doped with RE®+ are grown by magnetron sputtering technique as described in refe¬ 
rences.!^®^ Those composite layers contain silicon nanograins (Si-ng) and RE^+ ions dispersed in non-stoichiometric 
Silicon oxide (SiOa,) as schematically presented of figure]^ Those Si-ng act as sensitizers of RE^+ ions in their 
vicinity by a mechanism of energy transfer and therefore enhance their excitability by several order of magnitude. 
The effective absorption cross section of RE®+ sensitized by Si-ng is of the order of 10“^® cm^ against RE^+ di¬ 
rect excitation absorption cross section which is in the 10“^® to 10“^^ cm^ range. According to our experimental 
investigation^^! an exciton is photogenerated in Si-ng by pumping of SiO^, : Si-ng :RE^+ layer with an electroma¬ 
gnetic wave at 488 nm. After non-radiative deexcitations, the exciton energy reaches the corresponding energy 
gap of RE^+. Exciton can then vanish either by transfer of energy to the RE®+ ions by dipole-dipole interaction 
or by radiative or non-radiative recombination to Si-ng ground level. 

Hence those active layers have typical/particular photoluminescence properties as presented on figure]^ The 
characteristic emission band of radiative transition occurring from excited (Si*) to ground state (Si) in Si-ng is 
at 770nm. The characteristic emission bands of Er^+ and Nd^+ ions are respectively at 1532 nm and at 921, 1106 
and 1395 nm. 


3. TWO LOOPS ADE-EDTD METHOD 

The FDTD method is based on time and space discretization scheme of Maxwell equations proposed by 
Ye^i^ which allows to calculate the propagation of electromagnetic field (E,H) in time domain.E^ The ADE 
method consists in the use of extra terms such as current density J or polarization density P which are solutions 
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Figure 1. Optical amplifier waveguide based silicon rich silicon oxide doped with rare earth ions : (a)in longitudinal view 
with configuration of co-propagation of pump and signal and (b) in cross section view 




Figure 2. Schematic representation of the principle of optical indirect excitation of the rare earth ions via sensitization 
by silicon nanograin in silicon rich silicon oxide sample with silicon nanograins and doped with : (a) erbium ion (b) 
neodymium ion 


















Figure 3. Typical photoluminescence spectra of silicon rich silicon oxide sample containing silicon nanograins and doped 
with rare earth ions : (a) erbium ions (b) neodymium ions 


of differential equations with the aim to model some non linear optical behavior such as dispersive or gain 
mediaElEll The Maxwell equations describing fields E, H are rewritten in the form given in equation 


( 1 ) 


In equation!^ and are respectively the static dielectric permittivity and magnetic permeability of different 
considered materials, a is the electrical conductivity and p is a fictitious magnetic resistivity used for boundary 
conditions of the calculation box. The perflectly matched layers (PML) scheme developped by Berengei^ have 
been implemented as boundary conditions. Pjot = ^ P^- is the sum of all polarizations corresponding to each 
transition of absorbing/emitting centers (hereafter : silicon nanograins and rare earth ions). The use of one 
polarization density P^ per optical transition between levels i and j allows the description of the active layer 
dynamic permittivity e(a;) arising from the dipole moment densities induced by optical transitions in emitting 
centers. The FDTD space step A, time step Aj, Courant number Sc and calculation box size are reported in 
Table □ 


Table 1. FDTD algorithm parameters 


Parameter 

A 


Sc 

box length x 

box length y 

box length z 

Value 

50 nm 

10'"'' s 

1 

169A 

321A 

161A 


Neglecting the Rabi oscillation termp^ for a transition between levels i and j the polarization density Py is 
linked to the instantaneous electric field E(t), to the population difference AA^- = Ni — Nj, and to a constant 
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through the Lorentz type polarization density differential equatiorfi^ : 


+ Au:,, AA,,(t)E(t) (2) 

The time evolution of electronic levels population is modeled by classical rate equations system. Full detailed 
rate equations system of considered emitters will be given thereafter. As mentioned in the introduction, some 
excited levels characteristic lifetimes may be as long as a few msElE3 Due to the time step of ADE-FDTD 
classical method (figure]^), which are lower than 10“^^ s,li^the number of iteration should be as huge as 10^^ 
in order to reach the steady states of the levels populations. A conventional calculation is therefore impossible 


























in a reasonable time. Consequently a new two loops algorithm based on ADE-FDTD method (figure |^) has 
been developed that overcomes the multi-scale times issue of such a system and allows to describe the spatial 
distribution of the electromagnetic field and levels population in steady state in an active optical waveguide00 
This algorithm works with a drastic reduction of number of iterations to reach steady states values of fields and 
levels population from about 10^® to 10® iterations. 



Figure 4. Flowcharts showing (a) the classical ADE-FDTD algorithm and (b) the new hybrid-algorithm with the al¬ 
ternation of the short times loop calculating electromagnetic field and polarizations and the long times loop including 
calculation of levels populations. 


The description by the Lorentz oscillator polarization density (Equation 2) of an absorption or an emission 
process occurring during a transition i to j implies the equality of absorption and emission cross sections and 
only homogeneous broadening of the transition linewidth. At the resonance frequency, each Lorentz oscillator 
linewidth Awy is related to the absorption/emission cross section by the equation]^ 


Kij 1 

eoc Aujij 


(3) 


According to equation]^ high absorption cross sections aij (typically greater than 10“®^ cm^) lead to a small 
linewidth Auiij (about 10^® rad.s’^) which impose a large number of iterations to reach a steady state. In order 
to calibrate the proper absorption cross section while keeping the number of iterations as small as possible, we 
exploit the superposition property of the polarization densities by using a number Np of identical polarization 
densities with larger AuJij. Hence, equation ^ becomes equation Q. 
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At resonance frequency for ojij = 3.8 x 10^® rad.s“®, the two methods lead to identical cross sections : aij{Np = 
1) = aij{Np = 1000) as represented on figure]^ By this method, the cross section calibration is performed by 
choosing the appropriate value of the number of identical polarization densities Np and the linewidth Aw^. 

Silicon nanograins (Si-ng) are modeled as a two levels system (Figure where the ground and excited levels 
populations (respectively and Ngi*) are given by the rate equations (Equation and [^. Due to a low 
probability of multi-exciton generation in a single Si-ng,^ we assume the excitation of one single exciton by 
Si-ng, therefore the Si-ng population will correspond to the exciton population. 




















































Figure 5. Absorption cross section as a function of the pulsation with a transition at 3.8 x rad.s ^ for two number 
of polarization densities Np and two linewidths AuJij 
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Figure 6. Excitation mechanism of (a) erbium ions and (b) neodymium ions in silicon rich silicon oxide containing silicon 
nanograins and rare earth ions. Radiative transitions are represented by straight lines and nonradiative ones by sinusoidal 
lines. 
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The Si-ng physical parameters are reported in table 2 (for further detail see referenced ancd) . 
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Table 2. Parameters levels of silicon nanograins 


j ->■ i 

Lifetime (s) 

uJij {10^^ rad/s) 

rad/s) 

Np 

K (cm 

[Si-ng] (cm ®) 

Si* Si 

50.10"“ 

3.682 

1 

2756 

I(F* 

HP 


Erbium ions levels diagram is presented on figure |^. We consider two non-radiative transitions and one 
radiative transition from level 1 to level 0 at 1532 nm. Moreover an up-conversion process from level 1 to levels 
0 and 3 can be modeled. The time evolution of Er^+ levels populations is described by the rate equations [7| to 


taken from referenceai^lEHl and reported in Table ^ 


10 The concentration of erbium ions was taken equal to 10^° cm"®. All parameters of erbium ions transitions are 
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Table 3. Erbium ions parameters 


j ->■ i 

3->2 

2->l 

2->0 

l->0 

Lifetime (s) 

0.1 X 10"“ 

2.4 X 10"“ 

710 X 10"“ 

8.5 X 10"® 

ojij (10“® rad/s) 


1.23 



Aujij (10“® rad/s) 


0.15 



Np 


1 



Cup {cm'^.s~^) 




5x10"^'' 


Neodymium ions levels diagram is presented on hgure [^. We consider three non-radiative transitions and 
three radiative transitions and neglect up-conversion phenomena. The concentration of neodymium ions is taken 
equal to 10^° cm"®. The time evolution of Nd®"*" levels populations is described by the rate equations 11 
Parameters of neodymium ions transitions are taken fronJUEll and reported in Table 
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Table 4. Neodymium ion parameters 


j i 

4^3 

3^2 

3^1 

3^0 

2^1 

1^0 

Lifetime (s) 

230 X 

1000 X 10"'’ 

200 X 10"*’ 

250 X 10"'’ 

970 X lO"""" 

510 X lO”""" 

ujij (10^^ rad/s) 


1.34 

1.77 

1.99 



Aujij (10^^ rad/s) 


0.67 

0.18 

0.11 



Nj, 


1 

1 

1 




We aim at calculating levels populations in steady states in a faster and more convenient manner. Therefore, 
the preferred method would be to analytically calculate the steady state populations by setting = 0 in each 
equation. Consequently, this method is applied in case of neodymium levels equations system. In erbium rate 
equations system a finite-difference method was used in order to reach the steady states of population levels. 
This calculation was performed using a time step ten times lower than the shortest lifetime (0.1 /is) of the rate 
equations system. The calculation time is then no longer negligible but does not rise up significantly the global 
calculation time of the two loops ADE-FDTD method. 

4. RESULTS AND DISCUSSION 

Poynting vector (R = E x H) has been calculated from time Fourier transform of electromagnetic field (E, 
H) computed by ADE-FDTD method. Both signal and pump were injected as guided mode in the waveguide. 
We determine the main component Rz of the guided pump (488 nm) and signal (1064 nm) power density of R 
in transverse and longitudinal section view (in W.mm"^) on figures and respectively. A decrease of the pump 
■^pump g^jQjjg waveguide is observed, which is due to the absorption at 488 nm by the silicon nanograins. 
Moreover, the signal at 1064 nm does not appear absorbed over the 7fim length of the waveguide. 

The two loops ADE-FDTD algorithm allows to compute tridimensional population distributions in different 
states. For a transition of interest, N^igh and Nionj represent respectively the higher and lower levels of the 
transition. We present the relative population difference (RPD) distributions in the waveguide in a 

longitudinal section view along the propagation axis for the transition in Si-ng, and the 1532 nm transition of 
Er3+ (figure]^ and the 1064 nm transition of Nd3+ (figure [To|) for a pump power equal to 10^ mW.mm'^. 

In case of Si-ng (figures left, [lo] left), the decreasing RPD profile with waveguide length is consistent with 
the pump profile Rf“'”^’(z) decrease figure]^ due to the absorption of the pump field by the Si-ng. For erbium 
ions (figures 1^ right), the RPD of Nd3+ ions shows a decrease with propagation length in the waveguide. This 
decrease is characteristic of the pump strong absorption due to the presence of the nanograins as shown in our 
previous paper .1^ The RPD of Er3+ ions remains positive (inversion population is realized) over a length of 1.5 
fj.m, beyond which it becomes negative witnessing the threshold effect occurring with three-levels system. For 
neodymium ions (figures [To| right), the RPD of Nd3+ ions shows also a decrease with propagation length in the 
waveguide which is also due to pump strong absorption by nanograins. The RPD of Nd3+ ions remains positive 












































Figure 7. Transverse section views of Z-component of the time average Poynting vector of the pump (on the left) and 
of the signal (on the right) in the middle of the waveguide. The pump wavelength is Ap=488 nm and the signal As = 1064 
nm. Injected powers are respectively 1 W.mm~^ and 1 mW.mm"^ 
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Figure 8. Longitudinal section views of Z-component of the time average Poynting vector of the pump (on the left) and 
of the signal (on the right) in the middle of the waveguide. The pump wavelength is Ap=488 nm and the signal A3=1064 
nm. Injected powers are respectively 1 W.mm~^ and 1 mW.mm"^ 
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Figure 9. Relative population differences along the direction of propagation of Si-ng (on the left) and for the erbium ions 
(on the right) in SiOs, :Si-ng :Er®+ system for a pump power equal to 1000 mW.mm"^ 
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Figure 10. Relative population differences along the direction of propagation of Si-ng (on the left) and for the neodymium 
ions (on the right) in SiO^ :Si-ng :Nd®+ system for a pump power equal to 1000 mW.mm"^ 


(inversion population is realized) along the whole waveguide length witnessing the four-level system behavior, 
where the level 1 is depopulated quickly to the ground level leading to :S> Ni. 


This co-propagating pumping configuration allows to study locally the waveguides with various pumping 
values in one calculation. Indeed from population levels distribution Ni(x,y,z) we can calculated the gross gain 


per unit length (dB.cm'^) at the signal wavelength by equation 16 
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N/iig/i and Njoii, are respectively the higher and lower levels of the considered transition and Gabs and Gem are 
respectively absorption and emission cross sections. For erbium ions, we use '^tigh = Ni and N;o^=No and for 
neodymium ions ^uigh = N3 and Njou, = Ni (figure 6). Assuming equal emission and absorption cross sections, 
the gross gain per unit length longitudinal distribution for both rare earth ions are plotted in figure [m In case 
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Figure 11. Gross gain per unit length (dB.cm"^) along the direction of propagation in the case of erbium ions (on the 
left) and neodymium ions (on the right) for a pump power equal to 1000 mW.mm"^ 


of erbium ions, figure [TT^ shows a positive (in red) gain per unit length only in the first 1.5 /im of the waveguide. 
Elsewhere, the gain per unit length remains negative. For neodymium ions, the figure [T^ shows a positive gain 
per unit length that remains everywhere positive in the waveguide. On these two dimensional distributions of 


































the gross gain per unit length, the local gross gain per unit length can be extracted from anywhere in the active 
layer. We choose to take the local gross gain per unit length at the beginning of the waveguide in the center of 
the active layer (x = 4.5 /im and y = 8.55 /im) and in the first cell of the active layer of the waveguide where the 
pumping flux was recorded. We perform calculations with different pumping power densities at this extraction 
position and obtain the local gross gain per unit length plotted for different pumping power densities in figure 

m 



Pump power density (mW.mm'^) 

Figure 12. Local gross gain per unit length as a function of the pumping power density for a waveguide doped with 
Nd®+(open circle) and a waveguide doped with Er®+(open square) recorded at x = 4.5 /rm and y = 8.55 /rm (benter of 
the XY section of the active layer) and z = 0 (beginning of the waveguide). Losses found by Pirasteh et aP^ and by 
Navarro-Urrios et aP^ are reported respectively in blue and red dashed line. 


For Er®+ doped waveguide, above a threshold pump power of 1550 mW.mm"^, a positive gross gain is reached 
which increases up to 2 dB.cm'^ for the highest pump power density simulated. In order to estimate the net gain, 
we must account for the background losses such as those found experimentally by Navarro-Urrios et al.P^ They 
found on comparable waveguide losses of about 3 dB.cm'^ at 1532 nm, making it impossible to reach a positive 
net gain in this range of pump power density. In case of Nd®+ doped waveguide, we find that the optical gain 
remains positive over the whole power range and it increases up to 30 dB.cm'^ for the highest pump power of 
10^ mW.mm'^. Taking into account experimental background losses of 0.8 dB.cm'^ found by Pirasteh et aP^in 
a similar system, we can estimate a net gain per unit length. Figure shows that for a pump power above 130 
mW.mm'^ the losses (dashed line) can be compensated leading to a net gain. The gross gain per unit length 
obtained for Nd®+ remains higher that the one obtained for Er®+ whatever the pump power density range. The 
gross gain per unit length difference between Nd^+ and Er®+ doped waveguides increases also with pump power 
density. Difference in Si-ng/RE transfer efficiency linked to the levels dynamics as well as to the difference in 
absorption/emission cross section^ {ctet =6.10^^ cw? against aMd =1.10^® cm^) may explain this gross gain per 
unit length feature. The low gain and a short length of positive population inversion obtained by modeling of 
Er®+ doped based waveguide on the broad range of pump power make impossible the achievement of an optical 
amplifier with this configuration which would compete with commercially available systems. The gross gain per 
unit length reachable with Nd®+, about one order of magnitude larger than the one obtained with Er®+, could 
lead to a significant amplification. Nowadays, there is no commercially available comparable optical amplifier 
based on Nd^+ emission bands. However, Nd®+ doped Aluminum oxide channel waveguide amplifiers developed 
by Yang et aP^ show a maximal internal gain of 6 dB.cm’^ for a pump power of 45 mW. 

5. CONCLUSION 

A new two loops ADE-FDTD algorithm has been developed that allows to describe the spatial distribution of 
the electromagnetic field and steady state population levels in an active optical waveguide. The multi-scale times 
issue of such a system consisting in crippling calculation duration due to large difference between FDTD time 









step and population level lifetimes has been overcome by this algorithm. Moreover, we proposed a calibration 
method of the i —> j transition linewidth Aujij according to the experimental cross section, making a possible 
comparison between experimental and theoretical studies. We apply our algorithm to a strip loaded waveguide 
whose active layer is constituted of a silicon rich silicon oxide layer containing silicon nanograins and doped with 
either Er®+ or with Nd^+ ions. Our method permits to determine a three dimensional distribution of gross gain 
per unit length in the waveguide in steady state. We have demonstrated that, to obtain a net positive gain per 
unit length in SiOa, :Si-ng :RE^+ based waveguide, the neodymium ions are more suitable than the erbium ions. 
For the latter, the theoretical maximum gross gain per unit length of 2 dB.cm'^ at 1532 nm (lO"* mW.mm'^) does 
not compensate background losses experimentally estimated to 3 dB.cm"^. On the contrary, the use of neodymium 
ions leads to a gross gain per unit length of 30 dB.cm'^ at 1064 nm (lO"^ mW.mm'^). Moreover the background 
losses are compensated above a pump power threshold of 130 mW.mm'^. This theoretical demonstration of a large 
gross gain per unit length for a Nd®+ doped active layer may justify further experimental work in order to achieve 
Nd®+ doped silicon based waveguide optical amplifier or laser. Further studies may be performed exploring other 
concentrations of rare earth ions and Si-ng, other rare earth ions and other pumping configurations in order to 
investigate the possibility of achieving larger gain. This method may be applied successfully to describe steady 
states of other kind of emitters (quantum dots, quantum wells. Dyes...) and in other configuration (VCELs, 
down-converting layers...). 
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